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Abstract 

We analyse the proton Compton-scattering differential cross section for photon energies up to 
325 MeV using Chiral Effective Field Theory (xEFT) and extract new values for the electric and 
magnetic polarisabilities of the proton. Our EFT treatment builds in the key physics in two different 
regimes: photon energies u < m n ("low energy") and the higher energies where the A(1232) resonance 
plays a key role. The Compton amplitude is complete at N 4 LO, 0(e 2 <5 4 ), in the low-energy region, and 
at NLO, O(e 2 5 ), in the resonance region. Throughout, the Delta-pole graphs are dressed with 7rN 
loops and 7NA vertex corrections. A statistically consistent database of proton Compton experiments 
is used to constrain the free parameters in our amplitude: the Ml 7NA transition strength b\ (which 
is fixed in the resonance region) and the polarisabilities ctEi and /3mi (which are fixed from data below 
170 MeV). In order to obtain a reasonable fit we find it necessary to add the spin polarisability jmimi 
as a free parameter, even though it is, strictly speaking, predicted in xEFT at the order to which we 
work. We show that the fit is consistent with the Baldin sum rule, and then use that sum rule to con- 
strain olei+Pmi- In this way we obtain a^i = [10.65±0.35(stat)±0.2(Baldin)±0.3(theory)] x 10 -4 fm 3 , 
and Pmi = [3.15=F0.35(stat) ±0.2(Baldin)=F0.3(theory)] x 10 -4 fm 3 , with X 2 = H3.2 for 135 degrees 
of freedom. A detailed rationale for the theoretical uncertainties assigned to this result is provided. 
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I. INTRODUCTION 



Compton scattering from the proton, 7p — > 7p, at photon energies up to a few hundred 
MeV, has proven an excellent tool to study the subtle interplay of effective low-energy degrees of 
freedom of hadrons, and of the symmetries and dynamics which govern them; see recent reviews 
for details (H-EH. It reveals the chiral symmetry of QCD and the pattern of its breaking by 
probing pion-cloud effects, and the properties of the A(1232) as the lowest nucleonic excitation. 
Since it tests the (real and virtual) excitation spectrum of the target, probing the two-photon 
response of a proton also complements the information available in the one-photon response 
(e.g. form factors). 

Since the earliest experiments on the proton |I-15], a particular goal has been an extraction 



of the electric and magnetic polarisabilities, otEi and (3mi, which parametrise the stiffness of the 
proton against deformation in uniform, static electric and magnetic fields. At photon energies 
Wi a b below 50 MeV the Compton cross section is completely dominated by the Born terms, with 
the Thomson cross section as the low-energy limit. The polarisabilities manifest themselves as 
corrections to this point-like cross section which grow as wf^, but there is only a narrow range 
of energy before faster-varying terms contribute significantly. At ~ 150 MeV non-analytic 
structure associated with pion photoproduction appears, and around the same energy the tail 
of the Delta resonance also becomes apparent. In order to extract the polarisabilities, therefore, 
we need a reliable description of these effects. In this paper we analyze Compton scattering 
data in such a framework in order to extract a^i and (3mi in a model-independent way. Our 
tool is Chiral Effective Field Theory, xEFT, the low-energy version of QCD [l6|-[25|. In it, 



theoretical uncertainties can be estimated using its power counting, a systematic expansion 
scheme for physical observables. As the baryonic sector of Chiral Perturbation Theory (xPT) 
including nucleons and the A(1232), it orders all interactions consistent with the symmetries of 
QCD (and in particular the pattern of its chiral-symmetry breaking) by an expansion in powers 
of "low" momentum scales in units of a "high" scale A at which new degrees of freedom become 
relevant. Calculating both tree and loop diagrams to a given order in this small parameter 
results in an amplitude that is truncated at a certain, prescribed level of accuracy. 

Indeed, Compton scattering provided one of the early successes of xEFT applied to nucleons. 
Bernard et al. j26|, 2?} calculated the polarisabilities at leading one-loop order, without explicit 



A(1232) degrees of freedom, to be 

«LO = 10/3 LO = 10e2 9l = U Q 1Q -4 f 3 (U) 

in remarkably good agreement with both the extractions of these parameters extant at that 
time, and those obtained in subsequent experiments. Furthermore the cross sections obtained 
from that theory were in qualitative agreement with the then-available data below about 
150 MeV jU, 28]. Subsequently, a large number of experiments 29-4(| have measured un- 



polarised differential cross sections at photon energies between about 60 and 400 MeV; see 
the thorough discussions in the reviews [ll-0]. However, we shall argue in this article that the 
database remains quite sparse in important regions, in particular between 150 and 250 MeV. 
In addition, there are unresolved consistency issues between data sets from different experi- 



ments. Therefore, plans to measure Compton cross sections at MAMI and HI7S 
very welcome. 
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Results obtained for olei and (3 mi using ;\EFT may be compared to direct determinations 



from fully dynamical lattice QCD which appear imminent |44H48j . Such simulations, in turn 



take advantage of the fact that ;\EFT reliably predicts the strong dependence of the polaris- 
abilities on the pion mass, so that results can be extrapolated from unphysical quark masses. 
Thus, Compton scattering provides an excellent example of how %EFT serves as a bridge be- 
tween cross sections measured in experiments and the non-perturbative quark-gluon dynamics 
underlying the physics of hadrons. 

Apart from its relevance in Compton scattering, the magnetic polarisability 0mi also con- 
tributes the largest error on the two-photon-exchange contribution to the Lamb shift in muonic 



hydrogen [49H51 



And, while we focus here on the case of the proton, this study also provides in- 
put to extractions of neutron polarisabilities from data on elastic Compton scattering data from 
deuterium; see e.g. Refs (H, [Hj. Small proton-neutron differences stem from chiral-symmetry- 
breaking pion-nucleon interactions, probing details of QCD. For example, Walker-Loud et 
al. 53] recently found that the biggest uncertainty in theoretical determinations of the elec- 
tromagnetic proton-neutron mass difference now comes from the contribution of ~ 0.47 ± 0.47 
MeV from this effect pi. However, such effects can be probed with confidence only if both 



experimental and theoretical uncertainties are well under control. 

As already mentioned, there has been a significant increase in the data available on 7p 
scattering over the last twenty years. Over the same period the EFT description of this process 
has been refined in several ways. The pion-nucleon-loop amplitudes are now known to one order 



higher |54j-l56j . where two short- distance contributions to a^i and (3 mi enter the Compton 
amplitude. The polarisabilities are then no longer predicted, but can be determined from 
low-energy Compton data, and the fit quality indicates the extent to which %EFT correctly 
cap tures the energy- dependence of the Compton amplitudes. This has been done in Refs. 57, 



58 



In another line of attack, the A (1232) was included as an explicit degree of freedom, 
applying Lagrangians developed in Refs. j59-62[. This allows for a description which applies 
from zero photon energy into the Delta resonance region, and thus also for using data at 
intermediate energies for alternative extractions of the polarisabilities [i 



63 . In both these 



variants, baryons have traditionally been included as non-relativistic degrees of freedom in a 
version called "Heavy-Baryon Chiral Perturbation Theory" (HB^PT) [19|, |20|, |26j. A good 



qualitative reproduction of Compton data in this regime has also been obtained in a ;\EFT 
calculation which does not invoke the heavy-baryon expansion [64l-67|. and the agreement 
between the cross sections with and without the heavy-baryon expansion is very good — even at 
leading one- loop order — once the polarisabilities are fixed 68]| . In line with the power-counting 
philosophy, it is encouraging that the extractions or predictions of chei and 0mi have not varied 
dramatically from the original predictions, as demonstrated by the values advocated in a recent 
review pj. 

In this paper, we merge these three refinements: higher-orders in the chiral sector, explicit 
Delta degrees of freedom, and partially covariant formulations. In low-energy Compton scat- 
tering, one identifies the expansion coefficients 



p = UhL 

A ' 



Ma — M 



N 



A 



and — . 
A 



(1.2) 



where A 700 MeV is a typical large scale. Eq. ( 11.2)1 indicates both that different power 
countings involving choices about the relative size of P and e are possible, and that different 
kinematic regimes exist, depending on whether u ~ orw~ — M^. In order to account 
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for the hierarchy of mass scales m n < Ma — Mn, we follow Pascalutsa and Phillips [63] and 
identify an expansion parameter 5 ~ e while counting P ~ 5 2 . We differentiate between two 
photon regimes: at low energies, we count u ~ m n and obtain an amplitude which is complete 
at 0(e 2 S 4 ), where the leading (Thomson) amplitude is 0(e 2 ), while around the Delta resonance, 
we count u ~ 5 and the amplitude is complete at order 5 relative to leading order. 

We determine a^i and (3mi by fitting to the proton Compton-scattering database established 
in Ref. [lj. To that end, we constrain the parameters of the Ml 7NA coupling by the data 
around the Delta-resonance, and then extract polarisabilities from the data at lower energy. We 
find that an excellent fit can be obtained by including one contact operator which corresponds to 
a short-distance contribution to the spin-polarisability Jmimi — even though, strictly speaking, 
this effect only appears at one order beyond that to which our calculation is complete. In doing 
this we are departing from the strict EFT expansion, but the same strategy is necessary at 
lower order where, once the Delta is included, otEi and /3mi must be promoted (62]. 

This article is structured as follows. First, we define the theoretical ingredients of our 
calculation: the Lagrangians with both a covariant and nonrelativistic A(1232) in Sec. Ill Aj 
the Compton amplitudes and the definition of nucleon polarisabilities in Sec. IHBj the power 
counting in the different kinematic regimes and how to arrive at an amplitude which is valid 
across regimes in Sec. Ill C( explicit results for the 7NA vertex corrections in Sec. IHDI the 
prescription for the kinematically-correct position of the pion production threshold in Sec. Ill E( 
and constraints from sum rules and other processes in Sec. Ill Fl Sec.QIIlis devoted to the fit itself. 
After discussing our fitting strategy in Sec. IIII Al we present our results in Sec. IHI Bl followed 
by a discussion of details, alternative scenarios and convergence issues as well as, importantly, 
an estimate of the residual theoretical uncertainties in Sec. IIII CI After the Summary and 
Outlook of Sec. IIVI an Appendix gives references for the explicit forms of the components of 
the amplitudes. Preliminary versions of our findings appeared in Refs. [j], S 



Wi- 



ll. THEORETICAL INGREDIENTS 



A. Lagrangian 



Since the xEFT Lagrangian is well-established, see e.g. Refs. [2l|, |25]; we quote here only 
the terms relevant to Compton scattering at the order to which we work, so as to establish 
notation and conventions. The leading meson-meson terms are: 



4 2 ) =1 (d^ a )(d^ a ) + eA^ e 3ij TC{ d^TCj + 1 e 2 A^{v\ + n 2 2 ) 



1 22 
-mjK + . . . 



(HI) 



Here, ir a is the pion isotriplet with isospin index a = 1, 2, 3, f n the pion-decay constant, A^ 1 the 
photon field, and e = — |e| the electron charge. 

In heavy-baryon ;\PT, the relativistic nucleon Dirac spinor and isospinor ip is re-expressed 
in terms of a "heavy-baryon" spinor H that has the straight-line propagation of a nucleon of 
mass Mn and four-velocity v removed from its dynamics, so that fH = H. In the rest-frame 
of the nucleus its velocity is v M = g 0>1 and its spin is S^ 1 = (0, cr/2), where <7j are the Pauli spin 
matrices 19, 20]. At leading order, the pion and nucleon couple via the axial coupling g A : 



C ( ^ = H^wD + g A u-S)H 1 



{11.2) 
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where = <9 M — ieZA^ is the gauge- covariant derivative, Z = ~(1 + r 3 ) the nucleoli charge 
operator in isospin space, r a the Pauli matrices in isospin space and 

?V = — V(T a d^ a + ee a3 V a 7r b ^ + ...)• (H.3) 

Lorentz invariance generates the first three terms of the next-order Lagrangian: 

£ ^ =R] { 2]^ ( ( * ' ~ ° 2 ~ ^ ' °' V ' M} ) + 4Cl ^ ( X ~ 27^ 2 ) (IL4) 
+ (c 2 - (v ■ uf + c 3 u-u- -^[S", S»)eF^((l + ««) + (1 + «M)r 3 ) } # + . . . 

with the electromagnetic field tensor. The coefficients of the other five operators are low- 
energy constants (LECs): the isoscalar and isovector quantities and /c^ are the sum and 
difference of the proton and neutron anomalous magnetic moments, and c\ to c 3 contribute to 
7rN scattering. The latter in particular encode the effects of higher-energy baryon resonances, 
and so have different values in theories with and without an explicit Delta. We use the following 
values: 

Cl = _0.9, c 2 = 0.2, 5 3 = -1.6, (II.5) 

in units of GeV -1 . These numbers are obtained from fits in Ref. 24| to pion-nucleon and NN 
scattering in the theory without an explicit Delta, 

ci = -0.9±g;|, c 2 = 3.3 ± 0.2, c 3 = -4.7lJ;g, (II.6) 

by subtracting the Delta-pole contribution of ±4g 2 , NA /9(M / \ — Mn) from c 2 and c 3 respectively. 

The 0(e 2 P) Compton amplitude also includes a tree-level 1/M^ pi ece from £^ which has 
a fixed coefficient due to the low-energy consequences of Lorentz invariance, i.e. due to the 
1/Mjm expansion, and is required for the correct low-energy behavior of the forward scattering 
amplitude: 

4n = - # f ^j [S", S v ]eF^v°D v ({l + 2k«) + (1 + 2/^)r 3 )tf + h.c. (II.7) 

There are additional terms in that contribute to the 0(e 2 P 2 ) Born amplitude. We do not 
list them here, but they are discussed in the Appendix. In addition, the pion-pole graph enters 
which involves the anomalous part of d£l (^0123 — 

= ** w^ pct - (us) 

At 0(e 2 P 2 ), divergent loops need to be renormalised by operators that occur in C^. Those 
relevant for Compton scattering are: 



4n CT = 2vre 2 tf t [| ( ^« + 5 ^\r 3 )g, v - \(5a { °\ + S^\) + (tfag + ^£)r 3 



F^F v p E. 



(II.9) 
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They translate into additional, energy-independent contributions, Solei and 8/3mx, to the electric 
and magnetic polarisabilities which are isoscalar or isovector. This defines the 7rN sector of the 
calculation. 

For the A (1232) resonance, our calculation uses two variants: a manifestly covariant form, 
and its heavy-baryon reduction. In the latter, the field A^ of mass Ma is an iso-quadruplet 
heavy-baryon reduction of a Rarita-Schwinger field j60[. Since we work to NLO in the domain 
where the photon energy and the Delta-nucleon mass splitting A M = Ma — Mn are comparable, 
u ~ A M , we need selected terms from the Lagrangian up to third order fn\ : 



C 
C 

c 
c 



HB,(1) 
A — 




HB,(1) 




ttNA ~~ 


u 


HB,(2) 


iefti / 


7NA — 


M N V 


HB,(3) 


eb 2 


7NA — 


2M2 



(11.10) 

(11.11) 

(11.12) 

[iyV^( T 3 [jD ^/+])Aj + H.c.], 

(11.13) 

where a (fi) is the index of the isovector (vector) coupled to the isospinor (spinor), and the 
angled brackets indicate isospin traces. The choice of the parameters (j^a, bi, b 2 and Ma will 
be discussed in Sec. II III The term of the axial coupling g nNA is of order u, and the b\ interaction 
of order ecu. Both the b 2 and D\ terms give rise to couplings which are of order eu 2 . Terms 
generated by the heavy-baryon reduction are omitted in Eq. (III. 131) . Unfortunately, early works 



—if • D + A M ) A ap 

(H^d u 7r a Al + {A a v )^d v 'K a H +...) 



hUs-D F^v p + F^v p S • /)) Aj + 



4M| 



on the systematic heavy-baryon reductions of xEFT with an explicit A(1232) j59l-[6ll. 71 . 72 
were not consistent in their definition of 61; signs and factors of 2 come and go. The convention 

In Eq. (1II.13j) we have identified 



used here is in accord with that of Refs. 
Ref. 



71 



u m u% 73 



S 6 6 = -26 2 - 

These are the pertinent structures in the heavy-baryon Lagrangian. As we discuss in 
Sees. Ill C 21 and III D| we perform a fully Lorentz covariant calculation of the Delta-pole di- 
agrams. To do this, we employ the relativistic Lagrangian used in the (^-expansion fej. With 
the relativistic Rarita-Schwinger spinor, the free, 7rNA and 7NA pieces are 



63 



pp,(i) 



r PP,(2) 

7NA 



3e 



(11.14) 



2M N (M N + M A ) 



^{ig u F^ - g El5 Fnd^i-Kd ^g M F^g Elb F^ , (11.15) 



where = |[7 M ,7l, = lil^^T}, 7s = i7°7 1 7 2 7 3 , and F^ = \^ vpa F p(T . The (leading) 
terms involving the electromagnetic field that emerge upon heavy-baryon reduction of this La- 
grangian are equivalent to those in Eqs. (1II.12j) and flII.13j) . provided we make the identification 



g M = h{l + M A /M N )/3, g E = b 2 {l + M A /M N )/3. 



(11.16) 



In this presentation, we use the relativistic variant of the 7NA couplings, but quote results in 
terms of b\ and b 2 , using this conversion formula throughout, since these are more widely used. 
The s- and li-channel Delta-pole diagrams then reproduce the results given in the Appendices 
of Ref. [63j. However, the Lagrangians of Refs. [63j, |74j do not include a structure that reduces 
to that with coefficient D\ in Eq. ( 111.13ft ; this point will play a role in Sec. Ill Dt 
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B. Compton amplitudes and polarisabilities 

The standard decomposition for the proton Compton amplitude involves six non-relativistic 
invariants: 

T — Ax {e'* -e) + A 2 (?'* •£)(?■ k') + L4 3 • (e'* x e) + \A A a ■ (£' x £)(?'* ■ e) 

+iA 5 a • [(?'* x k) (e • I') - e X I') (?'* • k)} + L4 6 a • [(?'* x £') (e • k') - (e~*x k) (e'* ■ k)} , 

(11.17) 

where and e (k' and e*') are the three-momentum and polarisation vector of the incoming 

(outgoing) photon, and k (k') unit vectors pointing in the direction of k (k'). The A^s, i = 
1, . . . , 6, are scalar functions of photon energy u and scattering angle 9, and are real below the 
7rN threshold. 

All data is reported in either the centre-of-mass (cm) or laboratory (lab) frames. The 
corresponding differential cross sections are 

da 



|$| 2 |T| 2 , (11.18) 



cm, lab 

where the spin-summed and averaged, squared amplitude is: 



\jf = I|A 1 | 2 (l + cos 2 ^) + i|A 3 | 2 (3-cos 2 ^) 

+ sin 2 9 Re [AA* 3 A q + (A* 3 A A + 2A* 3 A 5 - A\A 2 ) cos 
+ sin 2 # [-\A 2 \ 2 sm 2 9 + -\A±\ 2 (l + cos 2 9) + \A 5 \ 2 (1 + 2 cos 2 9) + 3b4 
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|2 



2' 1 2' 

+2Re[A*(A 4 + 3A 5 )} cos 9 + 2Re[A*A 5 ] cos 2 9] . (11.19) 



The frame-dependent phase-space factor $ in Eq. ( III. 1811 is given by 

1 M N _ a 1 < b 



<& cm = - ^; <£> lab = - ^. (11.20) 



Following Ref. 56], we use Compton amplitudes A\-A§ in the Breit frame, in which the 
photon does not transfer energy to the proton and the amplitudes are either even or odd in the 
photon energy. The Breit frame variables u and 9 are then: 

u= 2MNU , costf = l-^. (11.21) 

JIM^t 2u 2 1 } 

where 

s — ii ! 

t = 2w tab w{ ab (cOS #l ab - 1), V = Jj^- = -(Wi a b + w( ab ) (11.22) 

s and t being the usual Mandelstam variables. See Refs. [H, 75] for more details. 

The most famous manifestation of chiral dynamics in Compton scattering occurs in the 
nucleon scalar polarisabilities, a^i and 0mi- These parametrise the global stiffness of the 
nucleon's internal degrees of freedom against displacement in an electric or magnetic field; 
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see our recent review [l[ for details. They affect a Compton-scattering reaction in which the 
photon makes El — » El and Ml — )■ Ml transitions, respectively. In consequence, the low- 
energy expansion of the Aj's can be used to define these electric and magnetic polarisabilities 
a E\ and /?mi- They occur as coefficients of the uj 2 terms in Ai and A 2 . For the Breit frame, 



such an w-expansion gives (keeping terms up to 0(1/M^)) [56l. 158 



Ai(u,\ 



Ao(cjJ 



A,(cj,e 



AAoj,9 



A 5 (u,9 



A 6 (u,t 



4M3 
e 2 u 
2M| 



+ 



4M3 



1 + K 



1 + COS I 



l Hi 



K 



cos 6) u 2 

+4n(a E i + M i cos6)cu 2 + C(w 4 ), 
{p) (2 + k^)u 2 cos6 - 4tc(3 m1 u 2 + C(w 4 ), 



((l + 2/t (p) )-(l W p) ) 2 cos#) +Af 



-AlXU^^ElEl + 1EIM2 + (JM1E2 + jMlMl) COS 8) + 0(u 5 



e 2 u 
'2M 2 



1 W p) ) 2 + Af -4vru; 3 (7 



MIMl 



1MIE2) + 0(UJ" 



e 2 ui 
2M 2 



[1 + K 



e 2 co 
'2M 2 



+ AI + 4tt^7mimi + 0(u b ) 

1 + K^) + Af + A^ lE lM2 + 0(lu 5 ) 



with k^ the proton anomalous magnetic moment, and 

e 2 g A a; 3 (cos# — 1) 



AI 



Al 



-A\ 



e 2 g A uj 3 



8n 2 f 2 (ml -ty 



(11.23) 



(11.24) 



^ 2 f 2 (m 2 n0 -t) ' 

the contributions from the pion-pole graph (see Fig. IH.lf ii) below). In Eq. (1II.23j) a^i and 
(3mi are supplemented by 4 spin- dependent dipole polarisabilities Jxiyv which describe photon 
transitions of definite multipolarities XI -h- YV , with V — I + {0; ±1}; the relation to Ragusa's 
spin polarisabilities 7i_4 [55( is given for example in Ref. [l|. 

The polarisabilities used in Eq . (III. 23|) can be generalised to define dynamical, i.e. energy- 
dependent, polarisabilities [62, 76]. (These should not be confused with the generalised polaris- 
abilities that can be accessed in virtual Compton scattering.) This is equivalent to a multipole 
decomposition of the structure (i.e. non-Born) part of the Compton amplitude. As in all mul- 
tiple decompositions, a complete set of dynamical multipole polarisabilities contains the same 
information as the Compton amplitudes, Ai, but examining predictions for these functions 
channel-by-channel can provide insight into how physical mechanisms may manifest themselves 
in data. For further discussion of the definition and usefulness of energy-dependent /dynamical 
polarisabilities, see Ref. [l[. Here, however, we work with the full amplitudes rather than a 
truncated multipole decomposition. 



C. Power counting regimes 

As discussed in the Introduction, three typical low-energy scales exist in Compton scattering 
with a dynamical A(1232): the pion mass m n pa 140 MeV as the typical chiral scale; the Delta- 
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nucleoli mass splitting A M = Re[MA] — Mn ~ 290 MeV; and the photon energy oj. Each 
provides a small, dimensionless expansion parameter when it is measured in units of a natural 
"high" scale A ^> A M , m n ,u at which xEFT with explicit A(1232) degrees of freedom can be 
expected to break down because new degrees of freedom enter. Typical values of A are the 
masses of the u and p as the next-lightest exchange mesons (about 700 MeV), so that 

P = ^ - 0.2, e = Ma 7 Mn - 0.4. (11.25) 
A A ' 

While a threefold expansion of all interactions is possible, it is more convenient to approximately 
identify some scales so that only one dimensionless parameter is left. In order to assign a unique 
index to each graph, we therefore follow the 5-expansion of Pascalutsa and Phillips [63]. It 
notices that m n <C — Mn and takes advantage of the numerical proximity P w e 2 in the 
real world to define 

i.e. numerically 5 ~ e ~ P?. The approach then separates the calculation into two different 
regimes, depending on the photon energy. In regime I, u) ^ « 140 MeV is "low", so 
that one counts u ~ m n ~ 5 2 A like a chiral scale. At higher energies (regime II, u ~ A M w 
300 MeV), one counts u; ~ A M ~ 5 A. 

This power counting accommodates the fact that Compton scattering changes qualitatively 
with increasing energy. In regime II, the photon carries enough energy to excite a A (1232) 
intermediate state whose large width and strong 7NA coupling makes it dominate the ampli- 
tudes. At lower energies, the Delta should play a less pronounced role. This will be confirmed 
in the list of all contributions in the following sub-sections. 



1. Regime I: u> ~ m n 

The Compton graphs that contribute up to the order to which we work can be separated into 
the three classes of Figs. III. II to [TT731 References for the actual form of the resulting amplitudes 
are given in Appendix |A] 

We first discuss the graphs separately for photon energies uj ~ m^. For contributions without 
an explicit A(1232), the only expansion scale in this regime is P = 5 2 , and the power counting 
is that of HBxPT, with only even powers of 5 contributing (63| . 

The first class consists of the "Born" or "non-structure" amplitudes of Fig. III. II and contains 
no loops. All graphs are calculated in the gauge v ■ A = 0, so direct 7N couplings do not 
occur at lowest order. Thus in regime I, the leading, O(e 2 S ~ e 2 P°), contribution is the 
"seagull" diagram Fig. IH.l( i). with a vertex from (111.41) ; this gives the Thomson amplitude. 
Corrections at 0(e 2 5 2 ~ e 2 P) are: minimal coupling to the proton charge and magnetic moment 
from (1II.4P in diagram (ii)(a); 1/M N seagulls from £® in diagram (ii)(b), and coupling 
to the anomaly via d£l flH-81) in graph (ii) (c) . At 0(e 2 5 A ~ e 2 P 2 ), corrections enter from a 
7NN vertex from £^ in diagram (iii)(a), and finally also from the short- distance contributions 
5a ei, 5 f3 mi to the scalar polarisabilities in flll.9p . diagram (iii)(b). 

The second class of contributions comprises the pion-loop diagrams of Fig. III. 21 At leading- 
one-loop order in this regime, 0(e 2 S 2 ~ e 2 P), the diagrams of Fig. III.2( i) contain only the 
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contribution with typical size 




U rsj A 




e 2 5° (LO) 


e 2 5° 


(ii) (a) \ / (b) X/ U° 


e 2 5 2 




(m) (a)\^ j (h)X \^y fJ 


e 2 <5 4 


e 2 5 2 



FIG. II. 1: (Colour online) Tree diagrams that contribute to Compton scattering in the e-v = gauge, 
ordered by the typical size of their contributions in the two regimes ui ~ ~ 5 2 and ui ~ A M ~ 5, 
respectively. The leading-order contribution in a particular regime is indicated by (LO). The vertices 



are from: £^ (no symbol), £^ (square), £^ (triangle), £^ (diamond), (disc). Permuted 

and crossed diagrams not shown. 

contribution with typical size u ~ to ~ A 



,(2) 



,(3) 



('1) 



.(4),o 



(i) (a) 




(c) 



(d) 




e^ 1 




e 2 5 4 



e 2 5 2 



FIG. II. 2: (Colour online) Pion-nucleon loop diagrams , ordered by their typical size for oj ~ ~ 5 2 
and uj ~ A M ~ (5, respectively. Notation as in Fig. III. 11 Permuted and crossed diagrams not shown. 



vertices of Eqs. flll.lj) and were first calculated by Bernard et al. [54|. At the next order, 
0(e 2 5 4 ~ e 2 P 2 ), one insertion of leads to two types of contribution, depicted in Fig. 111.2( h). 
Those of diagrams (a) to (k) involve 1/Mn terms which are required by relativity and therefore 
would already be present in a leading-order covariant calculation of the 7rN loop amplitudes, 
while those of graphs (1) to (r) involve in addition the LECs which are NLO in both the 
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contribution with typical size 




~ A 


(i) w 


e 2 5 3 


eH- 1 (LO) 


(ii) ( a ) y- 

— 1 — 


<{** (b)^X^ wV-n/ 
' f ' if ' 7 r 


e 2 5 3 


e^ 1 



FIG. II. 3: (Colour online) A(1232) and Air loop diagrams, ordered by their typical size for uj ~ 
m n ~ 5 2 and u ~ A M ~ 5, respectively. Notation as in Fig. III. 11 and double line: A(1232); shaded 
blob: 7rNA couplings including vertex corrections, as detailed in Sect. Ill Dl and Fig. III. 31 below. The 
given order for (i) of 0{e 2 5~ l ) for uj ~ A M applies only to the leading coupling at the vertices; other 
contributions will contribute at higher order. Permuted and crossed diagrams not shown. However, 
the s- and u-channel Delta-pole graphs in (i) occur at different orders in regime II, as detailed in the 
text. 



covariant and heavy-baryon approaches: c\ to C3, and kS & >. A calculation of all of these 
contributions at non-zero photon energy was first provided by McGovern [56]. (In Sec. Ill El 
we discuss the prescription which puts put the pion-production threshold at the kinematically 



correct position [56l . |62| . ) While individual loop diagrams at each order are divergent, in 
dimensional regularisation the full amplitude contains only two divergences at order 0(e 2 5 4 ~ 
e 2 P 2 ). Their structure must match the short- distance contributions to the polarisabilities of 
Eq. ( 1II.9P and Fig. III.l( iii)(b). since there are no other short-distance contributions to proton 
Compton scattering at this level, and the renormalisation-scale dependence in the LECs must 



cancel that of the loops |54|, |56J. These LECs encode contributions to the polarisabilities from 
mechanisms other than soft pion loops or Delta contributions, i.e. from short-distance effects. 
In contrast, the spin polarisabilities are still parameter-free predictions at 0(e 2 5 A ) in %EFT. 

Finally, two- and higher-loop diagrams and contributions from higher-order Lagrangians 
enter not earlier than 0(e 2 5 6 ~ e 2 P 3 ), and are thus suppressed. 

The diagrams in the third class contain a dynamical A(1232), as listed in Fig. 111.31 We 
exemplify the power-counting using the tree-graph (i), which includes the A (1232) propagator 

S { °\u^m n )oc 1 . (11.27) 

Since 5 2 ~ u ~ A M ~ 5 in regime I, the A(1232) propagator scales as A^, 1 ~ e _1 ~ 

and not as ~ 5~ 2 . Tree- level Delta effects are thus suppressed by one order in 5 relative 
to the leading-one-loop pion-nucleon diagrams. Here, we made the argument using the heavy- 
baryon propagator for the Delta, as it is more transparent, even though we actually choose to 
compute the pole graphs using a relativistic propagator, as described in Sec. Ill C 21 

A similar argument gives the power counting for the ttA loop diagrams of Fig. III.3( ii). In 
regime I, the loop introduces a factor of P 4 , while the pion propagators and vertices scale as 
e 2 P -2 , combining to e 2 P 2 ~ e 2 S 4 . The A(1232) propagator reduces this by only one power of 
5. Therefore, both the tree-level graph, (i), and the one-pion- Delta loop diagrams, (ii), scale 



oop 1 

as e 2 S 3 . These diagrams were first calculated by Hemmert et al. 6l|, [72| in the heavy-baryon 
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FIG. II. 4: (Colour online) The leading-order contribution to the A(1232) self-energy £a in regime II 
(u) ~ A M ). 



approximation, and we use those results in our calculation. The difference between results for 
7rA loops with the relativistic propagator, and those in the heavy-baryon approach is higher 
order in the expansion in 1/Mn ~ 1/A, and is numerically small j68(, so the inconsistency 
thereby introduced is insignificant. 

Corrections to these Delta-contributions from C^na a dd a P ower of P ~ 5 2 and thus enter 
only at 0(e 2 5 5 ). There are therefore no corrections to the amplitude from A(1232) effects at 
0(e 2 S 4: ) in this kinematic regime. 



2. Regime II: u ~ M A - M N 

In this regime, the A(1232) propagates close to its mass-shell, and its non-zero width is 
predominantly generated by 7rN loops. The details of the 5 counting for this case were derived 
in Ref. j63[ and are briefly repeated here. The leading effect comes from re-summing the leading 
one-Delta-reducible diagrams via a Dyson equation to obtain a dressed Delta propagator, so 
that Eq. ( III. 271) becomes 

4°V ~ A M ) cc 1 . (11.28) 

A M + i ImE A — u 

Since the kinetic energy flowing through the A(1232) propagator in this region is E ~ 5 and 
the 7rNA vertex of Eq. (III. lift scales as 8, the leading self-energy contribution, Fig. 111.44 scales 
as 8 3 . Written as a function of Mandelstam s, its imaginary part is: 

t rv / M ( 9ttna \ 2 (y^ + M N ) 2 -m^ , 3fl , M m v 

MM*)] =~{ mi ) k 6(k), (11.29) 

where 9(k) is the Heaviside function and k the centre-of-mass momentum of the pion in the 7rN 
system: 

' [ s _ ( Mn + m n ) 2 ] [s - (M N - ^Kf\ 
4s 



Following Ref. 63], we keep only the imaginary part of Sa and absorb the real part into Ma 
and into the Delta wave-function renormalisation. Note that we quote the relativistic result for 
the Delta self energy, for reasons to be explained. 

After it is dressed by the self-energy (III. 261) . the Delta propagator (III.28I) scales as S^(E ~ 
A M ) ~ 5~ 3 . This is different from the scaling of the propagator (III.27I) for u A AJ because 



as u becomes comparable to A M , there is a region where the difference u — A M scales as 
\co - A M \ ~ A5 3 . With A m 700 MeV, this is numerically the case for \cu - A M | rj 50 MeV. For 
photon energies this close to the A(1232)-nucleon mass-splitting, all terms in the Dyson series 
are of the same order, Ea ~ S 3 , so that the resummation of the Delta width is mandated by 
the power counting. This conforms with the reasoning that regime II should be limited to the 
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region about the resonance position, with its width approximately given by the width of the 
resonance. 

The fact that, relative to regime I, the photon energy is enhanced by one order to u> ~ 5, 
and the Delta propagator is enhanced by two orders, changes the relative importance of the 
contributions to Compton scattering in regime II. The leading order is 0(e 2 5~ 1 ), and does 
not consist of the Born graphs of Fig. III.lt Du ^ of the s-channel tree-graph, Fig. III.3( i). which 
contains the leading Ml coupling bi and the dressed A(1232) propagator obtained by solving 
the Dyson equation with the 0(<5 3 ) Delta self-energy due to N7r loops, Eq. flll.29j) . 

The propagator of the corresponding u-channel diagram is not similarly enhanced but counts 
as S~ l , as in regime I. Like the ttA loops of Fig. 111.3( h). therefore, its contribution scales as e 2 S 
in regime II, i.e. it is N 2 LO. 

Since we already use the relativistic self energy, NLO (O(e 2 5 )) contributions to the Compton 
amplitude come only from the 7NA vertex corrections, symbolically included in the shaded blob 
of Fig. 111.3( 1). These appear as both loop and counterterm effects, including one E2 coupling 
b 2 , and will be discussed in Sec. Ill Dl 

The only non-Delta contribution to order O(e 2 5 ) comes from the "seagull" graph of 
Fig. III. 1( a). which is now demoted to NLO in regime II. All other diagrams containing only 
nucleons and pions enter at even higher orders: the Born terms of Fig. IILl( ii) and iV7r-loops 
of Fig. HL2T i) at 0(e 2 5) (N 2 LO); and the Born terms of Fig. ITLlTiii) as well as the iVvr-loops 
in Fig. IIL2( ii) at 0(e 2 5 2 ) (N 3 LO). At these orders, additional interactions of the Delta with 
nucleons, pions and photons are found, together with two-loop contributions. 

So far our discussion of power counting has been based on HB^PT. A strict application 
would give the peak of the resonance at u = A M , with oj the Breit-frame photon energy. The 
correct position, a/s — Mn = A M , would then be restored perturbatively as an expansion in 
1/Mjsj. The cross section in this regime is dominated by the Delta resonance, and is very 
sensitive to its exact location. At a minimum, those corrections which ensure the correct 
real part of the resonance energy must be resumed. Doing this means that we avoid purely 
kinematic — and hence artificial — shifts of the polarisabilities from one order to the next. (See 
Sec. Ill El for the analogous procedure for the pion loops.) Note that, even with the correct 
real part of the resonance energy ensured, the full EFT amplitude up to NLO in the resonance 
region includes a number of l/M^-suppressed vertex effects. Furthermore, despite the fact that 
a/s — Mn <C Mn at the resonance, the strong dependence of Sa(s) on k (see Eq. ( III. 291) ) means 
that the relativistic and non-relativistic widths differ by a factor of 2 (for a fixed value of g nNA ). 

An efficient way of dealing with these issues in practice is to use a covariant expression for 
the Delta-pole graphs, based on the Lagrangian (1II.14j) [63]. This automatically takes care of 
both recoil and vertex corrections, by including effects which are higher-order in the power 
counting than the order to which we are working; except for the correct kinematics this cannot 
improve our agreement with data, but — as long as the power counting is valid — it should not 
hinder it either. In the same spirit, we also determine the resonance parameters Ma and g nNA 
from the Breit-Wigner peak and width, the latter via the relativistic formula. 

Since the kinematic shifts are irrelevant at these energies in the 7rA loops, we retain the 



heavy-baryon formulation for those [6ll . 1 72 
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3. Unified description of regimes I and II 

While separate sets of amplitudes for regimes I and II are possible, matching them smoothly 
in the transition region would be cumbersome. It is more convenient to work with one set of 
amplitudes which is applicable in both regimes simultaneously. We therefore choose to include 
all diagrams which contribute to N 4 LO (0(e 2 5 4 )) in regime I or to NLO (O(e 2 5 )) in regime II, 
even if they contribute at higher orders in the other regime. This is indeed the set of diagrams 
discussed above, Figs. III. ll to lII.3| with a relativistically resummed Delta propagator for the tree 
amplitudes of Fig. III.3( i). Moreover, each set is gauge and renormalisation-group invariant by 
itself. Therefore, adding them does not lead to conceptual problems, but neither does it improve 
the accuracy of the calculation. It is merely convenient to describe Compton scattering on the 
proton with one set of amplitudes from zero energy up to the Delta resonance. The accuracy 
in each regime is still given by the highest order at which all contributions are included. In 
regime I (oj ^ m n ), that order is 0(e 2 5 4 ) or 4 orders beyond LO, and the relative systematic 
uncertainty of the Compton cross section can thus a priori be assessed as <5 5 ~ [| . . . |] 5 ~ 1 
to 3%. In regime II (u ~ A M ), that is O(e 2 5 ), i.e. NLO, with a relative a priori error of 
5 2 « [§ . . . |] 2 « 15 to 25%. 

Since the amplitudes are therefore less accurate at high energies, we choose in the following 
to analyse in detail the data lying in regime I, and then check the consistency with experiments 
in the Delta-resonance region. 

We finally stress that the regimes of applicability of the different counting schemes are es- 
timates: transition from one regime to the next is not abrupt but gradual, and the above 
considerations provide only typical a priori error estimates. Convergence of the xEFT expan- 
sion order-by-order must be checked carefully so as to reliably assess theoretical uncertainties. 
Ultimately, comparison to data will help to decide whether the theory accurately describes 
Compton scattering. 



D. 7NA vertex corrections 

As alluded to in Sec. Ill C 2\ corrections which contribute to the vertices of the Delta-pole 
graph of Fig. III.3f i) are, in the low-energy regime, of higher order than that to which we 
work. However in the medium-energy region they are promoted to start at NLO; see Fig. III. 51 
the vertex corrections are 0(eu 2 ), whereas the dominant (61) vertex is O(eco). Thus pion- 
loop effects in the 7NA vertex contribute to the O(e 2 S ) part of the Delta-pole graph. They 
are therefore predicted to be more important in the vicinity of the resonance than even the 
basic 7rN-loop graphs depicted in Fig. 111.21 Since we use xEFT to fit Compton data in the 
resonance region — especially to constrain the Ml 7NA coupling — we include 7rN vertex loops 
in our calculation in order to be complete to NLO in this region. They have not been considered 
in recent Compton scattering calculations [2, |66, 67], although the importance of the imaginary 



part of Fig. |II.5f ii) (b) was recognised in the original 5-expansion study of Ref . 63 

Inclusion of these graphs is also necessary to satisfy Watson's theorem j77[, which requires 
that the pion photoproduction amplitudes at the peak of the resonance are purely imaginary. 
This implies that the real contribution from 77rN Born diagrams must be cancelled by Delta- 
pole graphs with a purely imaginary 7NA vertex that results from 7rN loops. When these 
elements are assembled into Compton diagrams the net effect is that, at resonance, there are 
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contribution with typical size 




UJ A 


(i) 


e,5 2 


e 8 1 


(ii) (a)\ (b)^-^ 


e 5 4 


e S 2 




e S 6 


e 8 3 



FIG. II. 5: (Colour online) Contributions to the 7NA coupling of Fig. III.3[ ordered by their typical 
size for 00 ~ m 7T ~ 5 and 00 ~ A M ~ 5, respectively. Notation as in Figs. ILI] and IL31 7NA vertices 
from C^ A , Eq. 15|) . at 0(e5 2 ) proportional to g M (i.e. 61; square), and at 0(e<5 3 ) proportional to 
g E (i.e. 62; triangle). For diagrams (ii) (b) and (iii), the order at u ~ A M refers to the imaginary part 
only; the real part is demoted. Permuted and crossed diagrams not shown. 




FIG. II. 6: (Colour online) Pairs of diagrams whose imaginary parts (with cuts indicated by dotted 
lines) are related by Watson's theorem at the Delta resonance. Notation as in Figs. III. II and 111.51 

cancellations between the isospin-3/2 pieces of of the 7rN loop graphs of Fig. III.2l and Delta-pole 
graphs with imaginary 7NA vertices from the loops of Fig. 111.5( h) and (iii). Diagrams which 
are connected in this way are shown in Fig. 111.61 

Since we are using heavy baryons for the calculation of the loops of Fig. 111.2} but the 
covariant version for the Delta-pole graph, we have a choice about whether or not to use the 
heavy-baryon expansion in computing these vertices; neither choice is fully consistent, although 
the inconsistencies are associated with higher-order effects in 5-counting. In what follows we 
use the covariant formulation to accord with other aspects of our treatment of the Delta, such 
as the 7rN-loop which gives the width (see Sec. Ill C 2j) . 

Pascalutsa and Vanderhaeghen [74j computed both 7rN-loop graphs of Fig. III. 51 at leading 
order in the covariant formulation of baryon %EFT. Since their calculation was reported at 
the A(1232) peak as a function of Q 2 , while we are interested in the vertices for real photons 
away from the A(1232) peak, for completeness we give results for this case here. (Note that in 
Eq. (47) of Ref. |74j a 7 is a misprint for /3 T ) 

Both diagrams contribute to both the structures whose coefficients are h\ and 62 at tree 
level. These corrections are labelled as (i = 1,2) for the first diagram, Fig. III.5( ii)(b). 
since it contains the photon coupling to the pion; and as (i = 1, 2) for the second diagram, 
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Fig. III.5( iii). since in it the photon couples to the nucleon: 

'1-2/ 













- 2C*NA 




— Cna 





4 N) = -&< N ) + 2C NA / dyy/ dx [x(l-x-y)W + (l-x-xy)M N ]WM; 2 , (11.31) 
Jo ./o 

where PF = -y/s, and 

.M 2 = x(l - y)M 2 + (1 - x)m\ -x(l-x- y)W 2 - ie, 
M 2 b = (1 - x - xy)M^ + xml - x(l — x — y)W 2 - ie, 

CNA = (W^ (IL32) 



The imaginary parts of both diagrams can be calculated in closed form [741 ]; we have also 
checked them against Watson's theorem for photoproduction. In the heavy-baryon limit, the 
expressions for graph (ii)(b) (in which the photon couples to the pion in the 7rN loop) agree 



with those obtained by Gellas et al. (71|, and for (hi) (in which the coupling is to the nucleon) 
with our own heavy-baryon calculation — with the following proviso. In the covariant theory 
graphs (ii)(b) and (iii) are treated as being of the same (leading) order. However in the heavy 
baryon formulation (iii) is higher order, suppressed by 1/Mn, as is most easily seen in the gauge 
e - v — in which the leading photon-nucleon coupling is the magnetic one. At the same order 
the anomalous magnetic moment enters, so that the 7N coupling is proportional to the total 
magnetic moment. Since we work to the order at which the anomalous magnetic moment enters 
for the direct pion loops, we do the same for the vertex loops, and multiply b± by the isovector 
Dirac magnetic moment (jy) = 1 + . Furthermore, even at this order in the heavy-baryon 
theory, there is no contribution from graph (iii) to the electric coupling 6 2 , the leading term in 
an expansion of &£ N ) from Eq. (IIOT]) being 0(1/M 2 ). 



Following Ref. [74J we subtract from the above results the real parts at W = Ma, so that 
the renormalised Lagrangian constants 61 and 62 are defined to give the real part of the dressed 
7NA couplings at the A(1232) pole. In what follows, when we discuss fitting b\ we mean this 
quantity. 

However, even after such a subtraction, both contributions to b\ (W) contain log(M>j) de- 
pendences. This corresponds to a residual logarithmic dependence on the renormalisation-scale 
in the heavy-baryon case [7l|. That could be absorbed (for (ii)(b)) by the HB Lagrangian LEC 
Di of Eq. flII.13p . but no corresponding LEC was included in the relativistic Lagrangian of 
Ref. [74]]. Since the expressions are finite, and we have no information on the LEC Di, we do 
not consider the additional effect from it. 

In fact, in the vicinity of the resonance, the variation with W of the real part of each 
loops graph is formally of higher order than the imaginary part (the argument parallels that 
of Ref. (63j when discussing the self-energy of the Delta). Moreover, the b^s also receive (real) 
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FIG. II. 7: (Colour online) The effects of including loop corrections to the 7NA vertex; in all graphs 
blue is the real part and red is the imaginary part (which vanishes below W = Mn + m w ). Dotted 
lines show no loop effects, solid the (renormalised) contribution of diagram (h)(b) and dashed, both 
(ii)(b) and (iii) (imaginary parts only, with the // v ) factor for bi). 



contributions from diagrams with the same topology as diagram (ii)(b), in which the baryon 
propagating in the loop is a A(1232). A full computation of the vertex functions to this order 
is however beyond the scope of this paper. And indeed, at NNLO in this kinematic domain the 
Delta propagator also receives substantial two-loop corrections [78j]. Therefore in what follows 
we discard the real part of diagram (iii), even when, working to fourth order, we retain the 
imaginary part for consistency with Watson's theorem. 
To summarise, therefore, the dressed vertices we use are 

bf\w) = k + i lm[b\ v) (W)) + Re[b^(W) - &?°(M A )] 

bf{W) = bf\w) + i /i (v) Im[&f°(W0], bf(W) = bf(W), (11.33) 

where superscripts (3) and (4) indicate use in calculations accuratein the low-energy regionto 
0(e 2 5 3 ) and 0(e 2 5 4 ) respectively. 

The results of this calculation are shown in Fig. III.7I The first two panels depict b\ and bi , 
with the dotted line giving the constant (tree-level) result for comparison. When the effect of 
diagram (ii)(b) is added the solid line is obtained (so these are b\ 3 \w)), and the dashed line 
is the result when the imaginary part of graph (iii) is included as well (with the ^jy> factor for 
bx, giving b^(W)). We see that the hierarchy of effects predicted by the power counting (after 
renormalisation) indeed holds: 

h > b { ? ] (11.34) 

b 2 ~ bf ] > 4 N) (H.35) 

(The inverted hierarchy for the loop contributions to 61 has a double origin: the contribution 
of (ii)(b) happens to be small, while the large value of ^y 1 makes the contribution from (iii) 
unnaturally large.) As for the relative size of b\ and 62, the ratio of multipoles E2/M1 is 

H = l m , where ^.gfWjl, (IL36) 

since the operator multiplying 6 2 contains an extra time derivative. The quantity 62 (W) ls 
displayed in the third panel, where the different scale on the y-axis makes it clear that the 
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effect of 62 on observables is higher order than the constant leading term in b\. It was from the 
ratio (1II.36j) . evaluated at W = Ma, that Ref. (74J obtained the result 62/&1 = —0.34. We also 
use that result to fix the ratio b^/bi. 

We emphasise that the results we obtained for the variation in the real part of bi(W) with 
W are dependent on our particular definition of Delta-pole piece. In particular, a redefinition 
of the A(1232) field would not change the value of these functions at W = Ma, but could affect 
their variation as we move away from the pole. Therefore the results shown in Fig. II 1. 71 only 
have meaning within the context of the specific £ 7 na and C^na that we have adopted here. 



E. Re-summing recoil corrections 

In the heavy-baryon version of one-nucleon %EFT, interactions and amplitudes are expanded 
in powers of 1/Mn- As with the correct position of the the A(1232) pole discussed at the end 
of Sec. Ill C 2\ that of the one-pion production threshold is thus only reached as a perturba- 
tive expansion with higher-order terms. This has consequences for Compton scattering since 
the opening of the photoproduction channel induces pronounced cusps in the amplitudes, as 
discussed in detail in Ref. [62]. Pion production starts in Compton-scattering experiments at 



photon energy m n + ~ 150 MeV in the lab frame. In contrast, the HB amplitude neglects 
recoil by assuming a nucleon which is at LO infinitely heavy, and thus puts the threshold at 
photon energy m n ~ 140 MeV. The 0(e 2 P) HB^PT amplitude has a cusp at this point, and 
that at 0(e 2 P 2 ) is singular there. The consequences of the cusp are seen even for u below m % , 
so getting it in the right place is mandatory for an accurate description of data around the pion 
threshold. 

Various schemes exist in the literature to correct for these purely kinematic effects by a 
resummation of higher-order terms that is motivated by a modified power-counting in the 



vicinity of a threshold (see Ref. [54| for an early discussion). All schemes agree within the 
order to which the xEFT calculation is performed, with each effectively trying to expand the 
amplitudes not about u> = but about the kinematically correct threshold position. As in 
Ref. [62| we replace the photon energy, u, with the quantity co s , defined via: 

uj s = \/s — Mn (11.37) 

This has the advantage that it implies changing the amplitude for u — > only starting at u 2 . 

We perform this substitution in all s and u channel (i.e. crossed) diagrams containing Ntt 
loops, Fig. 111.21 and take u to be the Breit frame photon energy. With this implementation, the 
Breit frame amplitudes are still crossing-symmetric. At 0(e 2 5 3 ) this is all that is done. Note 
that a slightly different way of obtaining the correct position for the irN threshold was adopted 
in Refs. 156 -I 



At 0(e 2 5 4 ) we need to ensure a finite amplitude as we cross the 7rN threshold. This requires 
that we do not double count the 0(e 2 P 2 ) effects which were already included in the 0(e 2 P) 
result via the shift (III. 37|) . Thus, if the straightforwardly computed HBxPT loop amplitude is 
(in an obvious notation) T^(u) + T^(u), at C(e 2 5 4 ) we compute: 

M + r<4> - ( W . - W )W , (n.38) 
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where the superscript (1) indicates that only the 0(u 2 /Mn) piece of the frame- dependent factor 
u s — u is retained. This yields a finite result for the amplitude as it crosses the cusp, and that 
cusp is now in the correct place regardless of the frame in which the Compton amplitude is 



computed [5J, [56 



In contradistinction, Ref. |62j performed the replacement in the cm photon energy and only 
for those s-channel diagrams which show an imaginary part above threshold, weakly violating 
crossing symmetry. The 0(e 2 S 4 ) diagrams of Fig. III. 2l were not considered there. 

We note again that — as long as we are not close to the specific locations of analytic structure 
in the complex plane, namely the Delta pole and the 7rN cut — neither the resummation of 
recoil corrections to place the 7rN threshold in the right place, nor our use of a fully- covariant 
amplitude for the Delta-pole graphs, add any new dynamics. In the heavy-baryon expansion, 
all these kinematic effects reduce in either regime to corrections which are higher than the 
orders with which we are concerned. 



F. Constraints from other processes and sum rules 

Strictly speaking, there are two parameters at 0{e 2 5 i ) which can only be determined from 
proton Compton scattering: the finite part of the short- distance contributions 5a ei, SPmi to 
the proton's scalar electric and magnetic dipole polarisabilities. The 7NA couplings 61 and 
62 can in principle be determined from the width and E2/M1 ratio of the A(1232) radiative 
decay. 

It is, however, useful to check the consistency of the values found for these parameters 
against other constraints. For example, the Baldin sum rule relates the sum of the scalar dipole 
polarisabilities at zero energy to an integral of the weighted total dissociation cross section 
ct 7 n^x starting at the production threshold uj t ^ r via a dispersion relation. We use the value 
from the recent global fit of (33| : 

am +Pmi = t^ [ d uj^^L = 13.8 ± 0.4, (11.39) 

where the canonical units, 10 -4 fm 3 , for scalar polarisabilities, are to be understood here and 
below. The "recommended value" of Ref. 0], 13.9 ± 0.3, which is based on an analysis of total 
photo-absorption cross sections, differs from the value ( III. 391) by less than the uncertainty in 
the sum-rule evaluation. 

Ref. [ll] also obtained a value for the backward spin-polarisability ■y 7T : 

7 W = -36.1 ± 2.1(stat) =f 0.4(sys) ± 0.8(mod), (11.40) 

in the units of 10 -4 fm 4 that are used throughout this work for the spin polarisabilities. Note 
that Eq. ( 111.40)1 includes the pion-pole contribution to 7^, in contrast to the way in which mul- 
tipole spin polarisabilities were defined in Eq. ( 111.23)1 . If also separately reports the statistical, 



systematic and model-dependence errors. Again, other values, compatible within error-bars, 
have been obtained by other authors: Ref. [31] gives a "recommended value" —38.7 ±1.8 from 
fixed-t dispersion relations, while 7^ = —35.9 ± 1.8 ± 3.2 was found in Ref. [2| from a fit to 
Compton data that used a more sophisticated fixed-t dispersion-relation approach. 
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Meanwhile, a recent reevaluation of the dispersion integral for 70 



79 



gives: 



7o = -Ieiei - Imimi - 1e\M2 - 1m\E2 = -0.90 ± 0.08(stat) ± O.ll(sys). (11.41' 
For further discussion, see Ref. 



III. RESULTS 



The data sets we shall use to fit our free parameters are discussed extensively in our recent 
review, Ref. Il|. For the reasons explained there, in the low-energy region we include data from 
Refs. @, @, MM, but not from Refs. 0, 0. The Baranov 150° data @, EI} were also 



excluded. Above the photo pro duction cusp, the data of Hallin [31| rises much more strongly 
than that of Olmos de Leon [33j . Hence the former was excluded from the analysis we presented 
in Ref. [l| . This leaves a gap in the database between co> lab = 164 MeV and wi a b = 198 MeV, so 
although we quote a cutoff of 170 MeV on the data set for our low-energy fits, in practice we 
could have chosen this to be anywhere in the range 164-198 MeV without altering the results. 
(For further discussion see Sec. IIII C 2l below.) Below our cutoff two further points are excluded 
as clear outliers, namely Federspiel (135°, 44 MeV) and Olmos de Leon (133°, 108 MeV). For 



the other Olmos de Leon points we follow Wissmann [80| and add a point-to-point systematic 
error of 5% in quadrature with the statistical error for the Olmos de Leon data set (not shown 
in plots). Lastly, although we remove the data of Hallin [3l| above 150 MeV, the data from 
this experiment at or below the photoproduction cusp are completely consistent with the rest 
of the database, and so, to maximise our statistical power, we include it in our fits. 

We float the normalisation of each of these data sets within the quoted normalisation un- 
certainty, by adding a piece to the usual x 2 function: 

2 2 _i_ 2 

"X. Xstat ' Xsyst. 

E' scts ^ / A/" i (dq-^/d^) eX pt - (dq^/df]) thcory \ 2 y^ s / A/} - 1 \ 2 , . 
, =1 h v a/}a.. ) u\ m J ' 1 j 

where (d<Jij/dQ) expt and Ay- are the value and statistical error of the ith observation from 
the jth experimental set, 8j is the fractional systematic error of set j, and Hj is an overall 
normalisation for set j. The Mj are additional parameters, which are fit by minimising the 



combined \ 2 . For more details, see Ref. 81] . (The MacGibbon data 32] is counted as two data 
sets since the higher-energy untagged data is treated separately). 

For the physical parameters in the 7rN sector of the theory we use 82] = 139.6 MeV, 



f n = 92.21 MeV M N = M p = 938.3 MeV, g A = 1.27, = -0.22 and «W = 3.71. The 
neutral pion mass m n o = 134.98 MeV is used in the tt° pole dia gram but not elsewhere; other 
isospin-breaking effects are neglected. The Delta parameters are fit to the Breit-Wigner peak 
and relativistic width as — Mn = 293 MeV and j TNA = 1.425. 

Finally, we remind the reader that we use units of 10~ 4 fm 3 for otE\ and (3 mi, and 10 -4 fm 4 
for the spin polarisabilities. 
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A. Fitting strategy 



The aim of this investigation is to extend and merge the approaches of Hildebrandt et al. 



and Beane et al. |57|, |58|, [62J, carrying out a definitive fit of olei an d 0mi to low-energy data at 
0(e 2 5 4 ) in chiral EFT with a dynamical Delta. As summarised in Sec. Ill C 3[ the ingredients 
are as follows. Born and pion-pole graphs are included, as are counterterms for Schei and S/3mi) 
see Fig. III. II In addition we add the 7rN and 7rA loops of Figs. III. 21 and 111.31 The Delta-pole 
graphs of Fig. III.3( i) are treated covariantly, using the Lagrangian of Eq. (1II.15j) together with 
the corresponding expression for the width Eq. fill. 29 j) . and we include vertex corrections as 
discussed in Sec. Ill Fig. 111.51 The amplitude is complete to (D(e 2 5 4 ) in the low-energy region 
(sometimes referred to below as "fourth order", since it corresponds to 0(e 2 P 2 = e 2 6 A ) for the 
7rN loops) and O(e 2 5 ) in the resonance region. 

Hildebrandt et al. additionally fit the 7NA coupling constant b\ to data from Refs. j3l|, [33[ 
which satisfied u < 240 MeV. This was done without the inclusion of any width for the A(1232). 
We have examined the advisability of repeating this strategy (with a Delta width included). 
In fits to data below about uj = 180 MeV, the 1-a-surface in the acEi-pMi-bi parameter space 
is large, with a particular degeneracy between b\ and cxei — &mx- Extending the fit to higher 
energies tends to constrain bi, but at the expense of preferring values of oei and (3 mi which 
are less compatible with the low-energy values. The significance of this finding is hard to judge 
in view of the fact that the data set between 170 and 240 MeV is sparse and inconsistent (as 
discussed in detail in the review [l|). 

Instead we adopted the following strategy. It clearly makes sense to use resonance-region 
data to determine 61; the cross section there is essentially proportional to bf. Unlike Hildebrandt 
et al., our amplitude is at least NLO in this region, since we use the covariant Delta. This 
further has the advantage that the parameters bi and 62 have the same meaning as in the 
photoproduction study of Pascalutsa and Vanderhaeghen [74j; we can therefore take the ratio 
= —0.34 from their work. Furthermore, though we choose to treat b\ as a fit parameter 
we should get a value which agrees with Ref. [lH {b\ = 3.81) up to corrections beyond NLO. 



However, above 200 MeV the data sets of, on the one hand, Hallin [31] and Blanpied 40], 



and, on the other, Wolf [3§J and Camen [39|, together with others working at Mainz |3j 
are too discrepant for a consistent fit to be obtained. We have chosen to use the Mainz data 
for our fits to the higher-energy region, for reasons discussed further in Ref. [2]. Thus we fit b\ 
to Mainz data between 200 and 325 MeV, then cxei and /3mi to world data below 200 MeV, 
and iterate until convergence is reached. We do not pay too much attention to the quality of 
the high-energy fit, since our amplitude is only NLO in this region; in addition the data are 
noisy and there are apparent trends in the data that are not fully reproduced. We do, however, 
expect a high quality of fit at low energies. 

In Ref. [ll] we presented an 0(e 2 5 3 ) fit to the Compton database using just this strategy (note 
though that 7NA vertex corrections were not included in the amplitude there). The Hallin data 
were ultimately excluded because of their unacceptable x 2 . We obtained a good fit with a x 2 
of 106.1 for 124 degrees of freedom (d.o.f.) (135 data points, 9 floating normalisations and 2 
fit parameters). This calculation with the 7rN-loops computed to third-order — 0(e 2 P) — in the 
xEFT expansion — sets a standard for the calculation we present here. 

However simply repeating such a fit with the 7rN-loop calculation improved to fourth order 
does not yield an acceptable x 2 . An examination of the work of Refs. 57|, |58|, |62J reveals why 
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this is the case: both the 0(e 2 P 2 ) pion loops and the Delta significantly improve agreement 
compared to the pure third-order HB%PT result by raising the cross section in the cusp region. 
Their combined effect raises the cross section above the data. The values of the spin polarisabil- 
ities also suggest this problem. They are significantly closer to the indicative values provided 
by dispersion-relation (DR) analyses in the Delta-less theory at fourth order in HBxPT than 
they are at third order. However, adding the Delta-pole then provides a substantial additional 
Imimi contribution j72j. This suggests a resolution of the problem: just as adding the Delta to 
the 0(e 2 P) Delta- less theory required the promotion of the fourth-order counterterms for ole\ 
and 0mi, so, when fourth-order 7rN loops and an explicit Delta are both considered, one or more 
of the (nominally 0(e 2 P 3 )) counterterms for the spin polarisabilites needs to be promoted, and 
included in our 0(e 2 5 A ) calculation. 

We have therefore explored the effect of varying each of the spin polarisabilities in turn, i.e. 
fitting one spin polarisability, along with ole\ and (3mi to the low-energy data. (Fits where 
more than one spin-polarisability is varied produce soft directions in parameter space which 
limit our ability to extract ole\ and (3 mi with reasonable errors.) We determine the optimal 
spin polarisability to include amongst the fit parameters according to the following selection 
criteria: the low-energy \ 2 should be good; the fitted value of the spin polarisability should be 
"sensible" — which in practice means bringing 70 closer to —1 than with the pure fourth-order 
fit; the Baldin-constrained fit should be compatible with that obtained without the constraint; 
and the reproduction of the high-energy data should not be badly compromised (for reasons 
already explained, we do not rely solely on the x 2 value to determine this). We find that on 
almost every count varying 7mimi is superior to varying ■Jeim2, 1mie2, or Jeiei] i n particular 
it gives the lowest low-energy x 2 ■ Furthermore Jmimi is the spin polarisability for which the 
pure fourth-order prediction is most discrepant with DR estimates 0], and the only one (with 
an acceptable x 2 ) with a fitted value closer to the DR estimates than the unfitted one. For 
all those reasons we choose to fit jmimi as an additional parameter when presenting our final 
0{e 2 5 4 ) result. 

B. Results with the 0(e 2 5 4 ) low -energy amplitude 

We fit both without and with the Baldin-sum-rule constraint ole\ + Ami = 13.8 ± 0.4 (see 
Sec. Ill Fj) . The statistical errors are obtained from the Xmm + 1 ellipsoid (ellipse). Since b\ is fit 
to an independent higher-energy data set, its small statistical error in principle feeds through to 
produce an additional small error on the low-energy parameters, but this is negligible compared 
to other uncertainties. 

Our result then is, for the fit without the Baldin sum rule constraint: 

ole\ = H.7 ± 0.7(stat) ± 0.6(theory), f3 M i = 3.8 ± 0.55(stat) ± 0.6(theory) (III.2) 

with 7mimi = 2.9 ± 0.65(stat) and b x = 3.62 ± 0.02(stat), as well as a x 2 of H0.5 for 134 d.o.f. 
For the Baldin-constrained fit, we obtain oei — 0mi — 7.5 ± 0.7(stat) ± 0.6(theory) or 

a E1 = 10.65 ± 0.35(stat) ± 0.2(Baldin) ± 0.3(theory), 

M1 = 3.15 =p 0.35(stat) ± 0.2(Baldin) =p 0.3(theory), ^ ' ' 

with 7mimi = 2.2 ± 0.5(stat) and 61 = 3.61 ± 0.02(stat), with a x 2 of 113.2 for 135 d.o.f. 



If we instead take b\ = 3.8 from Ref. [74|, the Baldin-constrained value of /3mi drops only 
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FIG. III. 1 : (Colour online) The projections of the la ellipsoid and Baldin-constrained ellipse onto the 
(Xei-Pmi plane, giving an ellipse and a line respectively. The grey band shows the Baldin constraint. 

slightly from 3.15 to 3.0, well within the statistical uncertainty. In both parameter sets we have 
also included estimates of the individual residual theoretical uncertainties from higher-order 
corrections. Before we can justify them in Sec. IIIIC44 we need to explore the sensitivity of the 
fits to choices of the parameters and data sets, as well as convergence issues. 

The statistical errors are calculated from the boundaries of the Xmin + 1 ellipsoid or ellipse. 
The projections of the la regions in parameter space onto the ole\-$m\ plane are shown in 
Fig. IIII.lt together with the Baldin-constraint band. We can see that the two extractions are 
fully compatible with one another, as we had required for an acceptable fit. We have also 
checked that the \ 2 is not unacceptably large for any individual data set, and that none of the 
normalisations floats beyond the quoted systematic error. In particular, the 65 Olmos de Leon 
points are fit with a \ 2 of 69 (unconstrained) and 71 (constrained) with negligible floating of the 
normalisation. Our slightly restricted set of world low-energy data is clearly highly compatible, 
and indeed the errors on some sets seem to be overstated. 



In Fig. IIII.2I the results are shown (along with third-order fits close to those of Ref. |l( which 
will be discussed later), together with all the world data, including those data sets that we did 
not include in our fits. As expected, the differences between constrained and unconstrained 
fits are small. Essentially no systematic deviation from the trend of the low-energy data is 
visible. The trend of the high-energy Mainz data is well captured, except for a tendency to fall 
somewhat low at forward angles. 

For reference, with b\ = 3.61 and 62 = — 0.34&1, the other three spin-polarisabilities take the 
values Jeiei = — 11, Iexmi = —0.4 and ^mxei = 1-9, (excluding the n° pole contribution which 
is -45.9 for 7^). See table 4.2 of the review |l| for values in other EFTs and in dispersion-relation 
extractions. 
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FIG. III. 2: (Colour online) Results of fits to Compton scattering data, as described in the text. The 
differential cross section in nb/sr is plotted as a function of lab energy in MeV, at fixed lab angle. Black 
and black dot-dashed curves are fourth-order Baldin-constrained and unconstrained fits respectively 
(with a band showing the 1 a errors on the former); green dashed (with band) and dot-dashed curves 
are the corresponding third-order fits (see Sec. IIII C 4p . Parameters are given in Eqs. (|III.3|) . (|III.2|) . 
(|III.5|) . (jIII.4j) respectively. The insets show the low-energy fit region, and the grey area is beyond 
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31, 
401. 



the high-energy fit region. The magenta/purple diamonds are MAMI data, principally Refs 
while the black squares are the data of Hallin & and the yellow stars are from Blandpied Ref. 
For the definition of other symbols see Ref. The data shown are within 5° of the nominal angle. 
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FIG. III. 3: (Colour online) Sensitivity of the fourth-order differential cross section (in nb/sr) to b\, 
plotted as a function of lab energy in MeV at three different fixed lab angles, 9^ = 45°, 90°, and 
155°. The solid line has b\ = 3.1, while the lighter band spans spans the range to b\ = 4.1. Other 
parameters are those of the fourth-order Baldin-constrained fit Eq. (|III.3|> . 



C. Further details of the fit 



1. Sensitivity of the differential cross section to b\, (Xe\, an d Pmi 

In Fig. IIII.3I we show the sensitivity of the low-energy fourth-order cross section to varying 
b\ within the rather wide band 3.1-4.1. Note that if 7a/imi is not fixed, it has a contribution 
from the Delta which is proportional to b\, and so the variation of the cross section with b\ is 
somewhat larger, especially at forward angles. However, the overall scale of the effect is similar. 
We observe that there is virtually no sensitivity to b\ below the cusp at more forward angles, 



which explains why Beane et al. |57|, |58[ were able to fit this region so well. 

In Fig. IIII .41 we show the sensitivity to varying olei ± (3mi within the bands 11.8-15.8 and 
6-9 respecttively. We see, as expected, that most sensitivity to chei + 0mi (&ei — Pmi) occurs 
at forward (backward) angles. Since there is no modern data below 60° it is not surprising 
that — as can be seen from Fig. IIII. 11 and as we will discuss further in Sec. IIII C 41 — sensitivity 
to + Pmi is significantly less than to chei — mi, suggesting that the Baldin-constrained fits 
are more reliable. 

Further exploration reveals that sensitivity to (3mi is almost vanishing at 90° even for energies 
up to the peak of the resonance. Of course the contribution to A% of Eq. (III. 23|) vanishes at 
this angle. However, there is a piece of 45 ~ \A 2 \ 2 , and this, in fact, is a dominant part at 90° 
in the resonance region. Nevertheless, the observation stands. 

From Figs. IIII. 31 and IIII. 41 together we can see the degeneracy between bi and — (3mi 
referred to in Sec. IIII A| and hence appreciate the difficulty of fitting both to the low-energy 
data alone. Our strategy of fitting b\ to data above 200 MeV eliminates this problem. 



2. Choice of data set 

After the elimination of inconsistent data sets, we are left with a choice of energy to define 
the top of the "low-energy" region. Accepting that Hallin data above the cusp cannot be 
accommodated means that there is no other modern data between 164 and 198 MeV, and 
indeed only Olmos de Leon data between the cusp and 164 MeV. Above 200 MeV not only is 
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FIG. III. 4: (Colour online) Sensitivity of the fourth-order differential cross section (in nb/sr) to ctEi± 
13 mi-, plotted as a function of lab energy in MeV at three different fixed lab angles, #i a b = 45°, 90°, and 
155°. The upper line shows the sensitivity to oei+Pmi', the solid line corresponds to cxei+Pmi = 11-8, 
and the band spans the range to 15.8. The lower line shows the sensitivity to ccei — Pmi] the solid line 
corresponds to oei — Pmi = 6, and the the band spans the range to cxei — Pmi = 9. Other parameters 
are those of the fourth-order Baldin-constrained fit Eq. (|III.3p . 



the data sparse and inconsistent, but the prominence of the Delta means that our amplitude, 
which is only NLO in the resonance region, is less trustworthy. So we prefer not to extend our 
"low-energy" fit into this region. On the other hand, since the amplitude above the cusp is 
markedly more different between third and fourth order than that below, a conservative choice 
of fit region would stop at 150 MeV — in which region the Hallin data is perfectly consistent. 
So to test our sensitivity to the choice of data sets and cut-off, we consider three possibilities: 
(I) Wmax = 150 MeV, Hallin included (137 points, 10 data sets); (II) w max = 164 MeV, Hallin 
excluded (135 points, 9 data sets); and (III) the union of both sets (147 points, 10 data sets). 

Fig. IIII. 51 shows that the results from all three strategies are highly consistent; the values 
for /3mi m the Baldin-constrained fits, with statistical errors, are (I) /3mi — 3.00 ± 0.36, (II) 
Pmi — 3.01 ± 0.36, and (III) (3mi = 3.15 ± 0.34 where excessive precision has been used simply 
to demonstrate the slight differences between the three. Though it clearly hardly matters which 
one we use, we choose as standard the largest data set, (III). 

We also consider the effect of changing the upper cut-off on the higher-energy data to which 
we fit bi; a cut-off of 350 MeV favours a somewhat larger b± of 3.69, and the Baldin-constrained 
/3mi m the corresponding low-energy fit drops to 3.08 (x 2 = 112), which is clearly well within 
our uncertainties. 
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FIG. III. 5: (Colour online) la curves in the PmvImimi plane for Baldin-constrained fits to three 
choices of data sets described in the text: Set (I) green, dotted; Set (II) blue, dashed; and Set (III) 
black, solid. (Ill) uses the most data and is our final choice. 



3. The effect of loop corrections to the 7NA vertices 

While the the 7rN loop vertex corrections of Fig. III. 51 have previously been considered in 
photoproduction [74]], this calculation is the first to include them consistently in a study of 
Compton scattering. Their principal effect is to generate a significant imaginary part in the 
vertex above the photoproduction threshold, and hence, since the Delta-pole contribution to 
the cross section scales as | t> x (^) | 2 , to decrease the resonance peak height for a given b\ (recall 
that Re(&i) is fixed to be the real part of the running coupling constant at resonance). In 
practice, this means that when we fit b\, we get a slighly higher value if the coupling runs than 
if it does not, somewhat reducing the slight discrepancy between the values obtained here and 



in Ref. [74j. A less significant effect is that the real part also runs, and we recall that this 
running is essentially arbitrary since there is a LEC (D\ of Eq. flII.13j) ) whose finite part we 
choose not to include. 

If we fit only to data below the cusp (where the imaginary part vanishes), the x 2 and Pmi in 
the Baldin-constrained fit differ by less that a percent with and without running. Including the 
rest of the Olmos de Leon data, the fit without running has x 2 — H7 and @mi — 3.2 ± 0.35 (cf. 
X 2 = 109 and /3mi — 3.15±0.35 with the running), which is only a very slight change. As shown 
in Fig. IIII.6[ once the low-energy parameters are adjusted, the curves below the cusp with and 
without the running are indistinguishable by eye. In the resonance region the effects are of 
course more marked, and there is a modest improvement in the fit when running is included. 



4- Convergence and residual theoretical uncertainties 

We are now in a position to examine the convergence of our results as we go from 0(e 2 S 3 ) 
'third order") to 0(e 2 5 A ) ("fourth order") in the low-energy amplitude. The nucleon Born 
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FIG. III. 6: (Colour online) Differential cross sections (nb/sr) as a function of cm. angle with (black, 
solid) and without (green, dashed) 7rN loop corrections to the 7NA vertex at three different energies. 
Best-fit Baldin-constrained parameters are used for each (see text). 



contribution changes appreciably between these two, whereas the difference between the fourth- 
order result and the fully covariant amplitude is negligible, so we choose to use the fourth-order 
Born for both amplitudes, just as we use the covariant Delta-pole graphs. (Since the fourth- 
order Born amplitude has polarisability-like terms, the main effect of this is to avoid a spurious 
shift in polarisabilities between orders.) Hence the only difference between the two orders 
(apart from the vertex corrections) is the graphs of Fig. |II.2f ii). These are of two kinds: 1/Mn 
corrections to the graphs of Fig. III.2( i). which would be considered part of the third-order 
calculation in a covariant formulation, and contributions from the second-order LECs q and 
which are fourth-order in both heavy-baryon and covariant formulations. Fig. IIII.7I 
shows the comparison between third order, fourth order with the LECs set to zero, and full 
fourth order, all with the same values of the polarisabilities a^i = 10.8, (3mi = 3. (Of course 
we do not switch off k in the Born terms in the first two calculations, only in the loops.) We do 
not show the effects of the q's and n separately, as the former is almost completely negligible. 
(Recall that in the theory with an explicit Delta the q's have to be readjusted, and are all 
rather small.) On the other hand k has an appreciable effect (largely due to k^'\ is small). 
The 1/Mn convergence is very satisfactory up to about 300 MeV, with the difference between 
the third and fourth orders only exceeding 10% where the third order is anomalously low, i.e. 
for forward angles at the cusp. We also see that the effect of k is, in the same range, of the same 
size as that of 1/Mn corrections. Hence there seems little to be gained in using a covariant 
formulation for the loops without including anomalous-magnetic-moment effects in the loop 
calculation too. Unsurprisingly, the situation is different in the resonance region, where the 
1/Mn corrections are becoming large. 

This convergence is very satisfactory from a formal perspective. However, quite a lot of 
the low-energy data lie in the region where the difference between the third- and fourth-order 
amplitudes is accidentally enhanced, and so identical results when one fits the polarisabilities 
to the data are not guaranteed. To test the convergence of the fourth-order results of section 
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FIG. III. 7: (Colour online) Differential cross sections (in nb/sr) at third order (green, dotted), fourth 
order without 2nd-order LECs (blue, dashed) and full fourth order (black, solid), all without 7rN loops 
in 7NA vertices and with a E \ = 10.8, Pmi = 3, b\ = 3.8 as well as unadjusted jmimi- Curves at 
w iab = 149 MeV (left panel) and wi a b = 260 and 340 MeV (right panel) are shown as a function of 
cm. angle. 



IIIIBI we perform a fit using the third-order amplitude, obtaining the following results: 

a E1 = 10.2 ± 0.5(stat) ± 0.8(theory), /3 M1 = 2.5 ± 0.55(stat) ± 0.8(theory), (III.4) 
with b\ = 3.69 ± 0.02(stat) and x 2 °f 118. 6 for 134 d.o.f., and for the Baldin-constrained fit 

CiEi = 10.7 ± 0.3(stat) ± 0.2(Baldin) ± 0.8(theory), 

p m = 3.1 =F 0.3(stat) ± 0.2(Baldin), ^0.8(theory), ^ ' ' 

with bt = 3.69 ± 0.02(stat) and x 2 of 120.3 for 135 d.o.f. 

The total x 2 is distinctly higher than for the fourth-order fit of Eqs. (lIII.3p . 2 j) but 
entirely acceptable; the preferred normalisation shift for the Hallin data is slightly greater than 
the systematic error of 4%, but omitting the data set has little effect on the fit parameters. This 
fit is shown together with the fourth-order one in Fig. IIII.2I In the review [l| we performed third- 
order fits without running couplings and excluding all Hallin data; for the Baldin-constrained 
fit identical values were obtained, while for the unconstrained fit the results were a = 10.5 ± 
0.5(stat) and (3 M1 = 2.7 ± 0.5(stat). 

The 1 cr curves for the third- and fourth-order fits together are shown in Fig. IIII.8I We see 
that the results are fully in agreement with one another. The range of ocei+Pmi accommodated 
in each fit is quite large, and they only just overlap at the l-a level. This can be explained by 
the fact, noted earlier, that the two cross sections differ most at forward angles near the cusp, 
which is also where the sensitivity to chei + 0mi is greatest, but there is little data at angles 

0lab < 60°. 

In contrast, the central values of ole\ — Pmi in both fits are extremely close, and neither 
changes appreciably when the Baldin constraint is imposed. In view of this, we prefer to quote 
the fourth-order Baldin-constrained values of Eq. (1III.3P as our final result. 
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FIG. III. 8: (Colour online) 1 a regions of parameter space for the third- and fourth-order, Baldin 
constrained and unconstrained fits. The green, dashed ellipse and the shorter error bar correspond 
to third order, the blue, solid ellipse and larger error bar are projections of fourth-order fits in which 
7M1M1 is also varied, which is why the regions are larger. The fourth-order regions are the same as in 

Fig. urn 



Finally, we can now return to the issue of the theoretical accuracy of the extractions at 
0(e 2 S 4 ) in Eqs. (IITL2|) and flITL3|) in Sec. fTTTBl Since the polarisabilities first enter at 0{e 2 5 2 ) 
with the parameter- free prediction of Eq. (II. II) . we expect corrections to be a priori of order 
5 3 ~ 6% relative to the LO result. With an average value (a^° + /3^)/2 m 7 to set the scale 
and 5 ~ 0.4, this leads to an error of ±0.5 (rounded up) for the individual polarisabilities of 
the free fit, or for the combination ole\ — Ami of the Baldin-constrained fit. 

To validate this naive dimensional estimate, we can check the convergence pattern from 
order to order, which is predicted to be 6 x 7 ~ 3 from LO to NLO and 5 2 x 7 ~ 1.1 from LO to 
NLO. For the free fit, a^i changes from LO to NLO by 2.3 units, and from NLO to N 2 LO by 
1.5 units, while (3mi changes in each case by 1.3 units, respectively. These are quite compatible 
with the power counting, and suggest an error on the fourth-order fit values of ±5 x 1.4 ~ ±0.6. 
To be conservative, this is what we quote on our non-Baldin-constrained fit flIII.2j) . 

However, Fig. IIH. 81 shows that the convergence for ccei — (3 mi seems to be much better than 
for oei ± 0mi, the latter being poorly constrained by the data. Judging by the shift from third 
to fourth order, the error on ctEi + Pmi could be taken to be as large as ±5 x 3 ~ ±1.2, while 
that on the oei — 0mi would be neglible. To estimate the error on the latter, therefore, we 
prefer to take S 2 times the LO-to-NLO change of 3.5, giving 0.6, and hence an error of ±0.3 
on the Baldin-constrained values of and ft mi separately. This also implies that the theory 
uncertainties for the free fit of ct^i and (3mi are highly correlated. 

A further way of estimating errors in an EFT is to look at the effect of the inclusion of 
partial higher-order terms. We have done this in the low energy region when we include the 
7NA vertex corrections of Section [II D[ The results for olei and (3mi with and without these 
corrections differ typically by about ±0.1. From Fig. IIH. 71 we might speculate that 7rN loops 
with two insertions of may in fact be the dominant correction at N 3 LO, but still this is a 
reassuring result which suggests that our theory errors are not underestimated. 

For the errors on the third-order results of Eqs. (1III.4P and ( jIII.5|) . we take the ±0.8 as 
estimated in the review [l[, values which are entirely compatible with the current analysis 
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though almost certainly too large for the Baldin-constrained results. 

We close by pointing out that the statistical and residual theoretical uncertainties in ex- 
tracting proton polarisabilities are now compatible. 



IV. SUMMARY AND OUTLOOK 



In this work we have presented a chiral effective theory analysis of 7p scattering data. We 
employed an EFT with explicit nucleon, pion, and Delta degrees of freedom, arranged according 
to an expansion in the small parameter 5 = m 7T /(M& — Mjy) = (Ma — M^)/A |63| . 

The xEFT amplitude accurate up to (D(e 2 5 A ) in the low-energy (u ~ m^) region was used 
(where the leading, Thomson, amplitude is O(e 2 S )). This includes the fourth-order 7rN loops 
of Ref. j56[, the it A loops first computed in Ref. j6l[, as well as Delta-pole graphs. At this order 
short-distance pieces of the proton scalar polarisabilities, ocei and Pmi, appear in the Compton 
amplitude and these, together with the 7NA Ml transition strength, b\, are the free parameters 
in our xEFT description of 7p scattering for u ~ m n . However, instead of fitting b\ in this 
low-energy region, we determine it using data up to u^h = 325 MeV and a xEFT amplitude 
which includes all O(e 2 5 ) (next-to-leading order) effects in the resonance region. Meanwhile, 
a E\ and {3 mi were fit to a statistically consistent Compton database that covers the domain 
< wi a b < 170 MeV. Fits to this database and to resonance- region data were iterated until 
convergence was achieved. This does not, however, result in an accurate description of data in 
the low-energy region. By including the short- distance piece of the proton spin polarisability, 
Imimi — nominally a higher-order effect in xEFT — in our 0(e 2 5 A ) amplitude we obtain an 
excellent fit to the low-energy data, jmimi is the only spin polarisability where such promotion 
yields a fit that is consistent with other constraints, and with the 0(e 2 S 3 ) fit presented in 
Ref. [lj . Since there is little 7p data at forward angles, and it is there that the differential cross 
section is most sensitive to (%ei + 0mi, we use the Baldin-sum-rule evaluation a^i + (3mi = 
(13.8 ± 0.4) x 10~ 4 fm 3 from Ref. 33j to constrain this combination of polarisabilities. 

The resulting fit has a x 2 of 113.2 for 135 d.o.f. We find, in the canonical units of 10 -4 fm 3 , 

a E i = 10.65 ± 0.35(stat) ± 0.2(Baldin) ± 0.3(theory), 

Pm = 3.15 =F 0.35(stat) ± 0.2(Baldin) =f 0.3(theory), (IV. 1) 

with 7mimi = [2-2 ± 0.5(stat)] x 10" 4 fm 4 and 61 = 3.61 ± 0.02(stat). The result for a E i - 0mi 
is very close to that extracted using the 0(e 2 S 3 ) amplitude, and is stable under modifications 
of the 7p database, and the inclusion of certain higher-order corrections to the 7NA vertex. 
Importantly, the value of b\ obtained in this fit is consistent with that found in the pion 



electroproduction study of Ref. |74| at the expected level of accuracy of both calculations 
(NLO). 

Eq. flIV.ll) presents a considerably larger value of ft mi than has been found in dispersion- 
relation extractions of proton polarisabilites from a similar database (see, e.g., Ref. (H). Since 
ctEi — ft mi is a free parameter in dispersion- relation fits, it would be interesting to know if 
such fits can tolerate the values of Pmi found here — and, if they cannot, what additional 
physics in the dispersion-relation calculation forces (3 mi to smaller values. We observe that a 
0(e 2 5 3 ) calculation in a variant of EFT that does not use the heavy-baryon expansion predicts 



Ami — (4.0±0.7) x 10 4 fm 3 66|, 67]. This is consistent with our extraction (llV.lj) . but disagrees 
with the dispersion-relation results in Refs. 0, [33[ • 
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The xEFT calculation in the resonance region could be improved to NNLO by including 
further corrections to the 7NA vertex and computing the pertinent pieces of the two-loop 
self-energy of the Delta in %EFT [78]. This would be particularly interesting for the develop- 
ment of a consistent xEFT description of pion-nucleon scattering, pion photproduction, and 
Compton scattering in the resonance region. The Baldin sum-rule value obtained in Ref. [33] 
should perhaps also be re-assessed, since slightly different values have been obtained by other 



groups [28l 1831] . The uncertainty on the value of a^i + (3 M1 obtained in this way is on the verge 
of becoming a significant component in the overall uncertainty in a^i and /3mi- Meanwhile, if 
the amplitude in the low-energy region were to be improved, a calculation beyond (9(e 2 5 4 ) in 
the 5 expansion would have to be carried out, and we note that graphs with two 7rN loops enter 
at 0(e 2 S 6 ). As the result f llV.lj) makes clear, at this stage the theory error is slightly smaller 
than the statistical error, and so it seems that the accuracy with which and (3mi are known 
can be more profitably improved by focusing effort in other areas. 

In particular, if more data were available between 7rN threshold and the A(1232) resonance 
then we could extend the low-energy fit region and obtain tighter constraints on oei and (3 mi- 
That data that exists here are sparse, and different experiments are not consistent with one 
another. In consequence our database has a sizable gap between wi a b = 164 MeV and u\ a h = 198 
MeV. Data of high precision in this kinematic domain, with well-documented systematic un- 
certainties, could help reduce the statistical errors on chei — 0mi- Indeed, if such experiments 
extended into the resonance region they might help to discriminate between the conflictin g re - 



sults of, on the one hand, Hallin l3l| and Blanpied |40j, and, on the other, Wolf (38J, Camen 
and other groups at Mainz 3414371] . 

Ongoing experiments at MAMI (Mainz) and HI7S (TUNL) to measure 7p asymmetries at 
u = 250-350 MeV (MAMI) and u « 100 MeV (HI7S), and thereby extract spin polarisabil- 



ities, will provide additional testing grounds for our approach |41|, |42j. Meanwhile, new data 
on absolute cross sections from 120-180 MeV will supplement the current database, and, if 
sufficiently precise, aid the extraction of o^i — [5 mi |43J. 

In fact, the xEFT amplitude developed here is highly successful in reproducing extant Comp- 
ton differential cross section data. Perusal of Fig. IIII.2I shows that it describes experiment over 
a wide range of energies and angles. This is possible because the amplitude contains the key 
non-analyticities in 7p scattering below oj\^ = 350 MeV: the pion photoproduction cusp and 
the A (1232) resonance pole. This benefit comes at a price, though: the chiral dynamics en- 
coded in these non-analyticities describes data so well that effects due to short-distance physics 
play a relatively small role. Consequently, we had to exercise significant care in the treatment 
of (sometimes discrepant) data, and also consistently incorporate a number of subtle effects in 
the 7p amplitude to obtain accurate values for ocei and (3mi- 
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Appendix A: Details of the amplitudes 

The nucleon Born diagrams calculated to fourth order are as given by McGovern [56j and 
reproduced in Eq. (III. 23[) : they are indistinguishable from the full relativistic form given by 
Babusci et al. [75] for the energies of interest. They are obtained from the diagrams of 
Fig. III.1( the third-order vertex represented by a triangle employs terms given in Eqs. (3.8) 
and (3.9) of Ref. the divergences represented by <i 24 and d 3 i in the latter renormalise 
in Fig. 111.2( h) (k) and the loop in (m). 

The fourth-order vertex represented by a diamond in Fig. IH.l( iii) (b) contributes to both the 
0(e 2 P 2 ) Born amplitude and to the polarisabilities. The LECs 8aEi and 5(3mi of Eq. fill. 91) 
result from linear combinations of the fourth-order operators Og$ to O94, Om and 118 of 
Ref. |§il|. The Born amplitudes also involve fixed-coefficient parts of these (table 6) together 
with those of the operators X 41 , Iw and Y n . The relevant /3-functions are those of Oi 6 to O 2 o, 
Oi52, 0i53, 0160 an d Oi98 of Ref. |85j (though there would seem to be a misprint in one or more 
of the "eye-graph" contributions as the divergence cancellation for SctEi does not quite work.) 

The pion-pole contributions are given by Bernard et al. j2l[ and in Eq. ( 1II.24|) ; this form does 



not change at fourth order and is indeed relativistically invariant, although at higher orders 
still, form factors will enter at the vertices and the LEC di$ shifts g A j ' f n to the numerically 
almost identical g xNN /M N . 

The 7rN loops of Fig. III. 21 are given at third order by Bernard et al. [2l| and at fourth order 
by McGovern [56] . At both third and fourth order, we use u s = y/s — Mn to shift the threshold 
to the correct energy, as detailed in Sec. Ill El The it A loops are given by Hildebrandt et al. [62j 
in their Appendix B, except that we use the Breit-frame photon energy u throughout in place 
of u s and u u to preserve crossing symmetry. 

The s- and ^-channel Delta-pole diagrams are calculated using the Lagrangian of Pascalutsa 
and Phillips [63}. The expressions for the Delta-pole amplitudes given in Eqs (A3) and (A4) 
of that paper refer to a redundant, covariant set of eight operators. We note that the correct 
forms of three of these are 

Ot = q»q' al av + l m qJ\ T = <f^ V + 1^'J\ <T = ^q a q' ; (A.l) 



the others are given in their Eq. (51). (Other misprints in Ref. |63| are noted in Ref. |l|; in 
particular the expression for the width given in their Eq. (42) is a misprint for our Eq. flll.29p .) 
These eight operators may be written in terms of the six Breit frame operators of Eq. (III.17P , 
which we will denote tj (with e.g. ti — e'* ■ e) using 

6 

e';ore u = 2M N ^ Qjtj, (A.2) 
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where 
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£ is the Breit-frame nucleoli energy, E = a/M^ — t/(4), and z = cos^ = t/(2w 2 ) + 1. 
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